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1.0 INTRODUCTION 


The tables in this report give the response of various types of capacitors 
to various transient gemma radiation environments. The tables were generated on 
the Burroughs B5500 digital computer at Stanford University. The parameters used 
to characterize the capacitors were obtained from the TREE Handbook (Battelle Me- 
morial Institute, DASA 1420, Feb., 1964) and "Analytical Methods for Predicting 
Transient Nuclear Radiation Effects on Electronics Circuits and Devices," (Boeing 
Corp., RTD TDR-63-3007, July 1963). ‘The capacitor parameters in these reports 
represent only the average response of many similar capacitors to many different 
rediation environments. Unit-to-unit and manufacturer-to-manufacturer variations 
in the capacitors make it impossible to predict the effects of radiation exactly, 
so these tables should be considered as giving only a "typical" response of the 
particular capacitor to the given radiation environment. Assuming, however, that 
the capacitor parameters from the above reports are exact, then it is believed 
that the numbers given in the tables, with a few exceptions discussed later, are 
accurate to three significant figures. 


Additional tables are given at the end of this report. These tables give 
the voltage left across a capacitor at various times after t=0 for a fixed resis- 
“tance shorted across the capacitor at time t=0. The purpose of these tables is to 
help the design engineer simulate the effects of radiation. The organization and 
use of these tables are described in later sections of this report. 


2.0 ORGANIZATION OF TABLES 
2.1 Capacitor Tables 


The one hundred forty-four capacitor tables in this report describe the re- 
sponse of six different types of capacitors to twenty-four different radiation en- 
vironments. The column headings in each table give the particular external circuit 
time constant for which the column was calculated, and the row labels at the left in 
each table give the time after the initiation of the radiation event for which the 
row was calculated. The numbers given in the tables are the voltages left on the 
capacitor under the conditions described in the row, column, and table headings, 
assuming the initial voltage to be one volt. Thus the organization of the capacitor 
tables may be depicted schematically as in Figure 1. 
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The radiation environments used are rather idealized versions of typical 
flash x-ray, nuclear veapon, pulsed reactor, or (with reservations about the 
effects of particle charge) pulsed accelerator transient radiation environ- 
ments. The general form of the assumed e'vironment is a short, high intensity 
rectangular radiation pulse followed by lower level steady radiation (see 
Figure 2). It is believed that some sort of approximation to the environ- 
ments listed above can be found either in the tables or by interpolating be- 
tween tables. The twenty-four radiation environments are found by taking each 
of the initial rates of 10'°, 1021, 10°”) ana 1025 Reds/sec ani for each initial 
rete using each of the final rates of 0, 10, 10°, 10°, 10°, and 10” Rads/sec. 
The 4 initial rates together with the 6 final rates give the 2h radiation en- 
vironments for which the capacitor tables are calculated. 


The initial pulse is assumed to last for 25 nanoseconds. This duration 
was chosen because: 


(1) 4t 1s convenient for numerical evaluations of the capacitor 
discharge equations, and 

(2) 41t approximately corresponds to the pulse length expected 
from the Lockheed Missiles and Space Company 2 Mev flash x- 
ray machine scheduled to be operational in April 1965. 


Normalization of the capacitor tables to different initial pulse durations is 
discussed later. 


The six types of capacitors chosen for inclusion in this report are those 
on which adequate information is available in the TREE Handbook. They are, in 
the order in which they appear in this report: 


Ceremic capacitors (barium titanate) 
Tantalua capacitors (solid) 

Mica capacitors 

Mylar capacitors 

Paper capacitors (Vitemin Q) 
041-Filled Paper capacitors (Guieman) 


¢ It is hoped that this report will be kept updated as information on other 
types of capacitors becomes available in the literature. In addition, present 

‘ plans call for characterization of certain types of capacitors, particularly 
glass and wet tantalum, on the IMBC 2 Mev flesh x-ray facility mentioned above. 


In almost all practical applications of capacitors, the capacitor will be 
receiving charge from an outside circuit while the radiation field is discharg- 
ing it. For this reason, the voltage left on the capacitor is calculated in 
these tables for a number of time constants of the external recharging circuit. 
These time constants are infinity (open recharging circuit), R.C = 10 sec, 1 sec, 
100 ms, 1( ms, 1 ms, 100 us, 10 us, and 1 us. These musbers form the column head- 
ings in the capacitor tables. In order to keep the tables to a reasonable, or in- 
deed, finite, size, it was necessary to make the following assumptions. 


(1) Figure 3 shows the equivalent circuit of the capacitor. ‘The 
switch is assumed to close at time t = 0. 

(2) The battery V, and the resistance R, represent the Thevenin 
equivalent circuit of the capacitor driving circuit during 
the radiation pulse, and do not change during the pulse; this 
circuit includes any inherent capacitor leakage resistance. 

(3) The battery voltage V, is the seme before and during the ‘Tedie- 
tion pulse, and the capacitor is charged to Ne before the pulse 
starts. 

(4) The resistance R, may be different before and during the pulse, 
but it must be possible to reasonably approximate it by a constant 
resistance during the pulse. The value of R, during the pulse must 
be used to calculate the time constant. 

(5) Permanent radiation damage is neglected. 


It is certainly true that most circuits will not exactly satisfy these as- 
sumptions during the pulse, but it is also true that for most cases these assump- 
tions are quite reasonable. 
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O 25 ns. Tins —> infinity 


Figure’2s Rediation pulse assumed in calculation 
of capacitor tables. 


Figure 3: poteciien eateries strexss The switch 
Gloses at the start of the radiation pulse. 
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The row labels in the capacitor tables are the times after the initiation 
of the radiation event for which the mubers in the row are calculated. The 
firet row is at time t = 0 for which the capacitor voltage is defined to be 1 
volt; the second row is at 25 ns, the end of the first pulse, and the third 
| ‘through tventy-fourth rows are times during the second rediation rate, going 
from 1 us to 10 sec in approximately logarithmic increments. 
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2.2 Resistance Simulation Tables 


The last nine tables are for the convenience of the designer in simulating 
the effects of radiation. 


3.0 USE OF TABLES 


3.1 Capacitor Tables 


Consider the following problem: Find the best capacitor for which date is 
available for use in a radiation environment consisting of an initial pulse of 10/4 
Rads/sec for 25 nsec followed by a constant rate of 10° Rads/sec. The time constant 


of the external charging circuitry is 10 msec and the circuit sust operate 100 msec 
, efter the start of the radiation transient. 


To solve this problem it is necessary to look up one data point for each type 
of capacitor. First look in that table for ceramic capacitors which is calculated 
for an initial pulse of 10 ana a constant rate of 10° Rede/sec. Looking at the 
intersection of the 10 msec time constant column witn the 100 msec time rov, we 
find that 99.2% of the initial capacitor voltage remains at 100 ms. Repeating this 
process for the other capacitors we get the following table: 


Capacitor type Final voltege/initial voltage 
Ceramic 0.992 
Tantalua 0.932 
Mica 0.997 
Mylar 0.852 
Paper 0.969 
01l-paper 0.826 


ae 


Thus ve would not want to use oil-filled paper or mylar in this application. 
Ceramic, mica, or possibly peper would be the indicated dielectric for this 
application, depending on the capacitance value needed and other circuit and 


permanent demage consideretions. 


As e second example consider the problem of interpolation. We might wish 
to interpolate between any combination of initial rate, final rate, time constant, 
or time. The two examples hgre will illustrate one possible technique of inter- i 
polation for one variable and for two variables. The technique illustrated may 
be used for any variable or combination of variables, but the amount of work in- 
volved goes up exponentially with the number of variables. 


As a first exemple, suppose we wish to find the voltage on a ceramic capaci- 
tor if the initial radiation rate is 107 Rads/sec for 25 nsec., the final rate 
is 10° Reds/sec, the time of interest is 100 msec after the start of the radiation, 
and the circuit time constant is 30 msec. There is a table for ceramic capacitors 
in this rediation pulse, and the voltages at 100 msec are tabulated. However the 
time constant of 30 msec is not included in this table. Thus it is necessary to 
interpolate to get the voltage for a 30 ms time constant. At 100 ms after the 
start of rediation the voltages left on the capacitor for various time constants 


Time constant voltage 
1 msec 0.993 

10 msec 0.932 

100 msec 0.505 

1 cec 0.141 


These mmbers are pictted (on a semilog scale) in Figure 4. If a smooth curve is 
drawn through the points, we find that the voltage left at 100 msec for a circuit 
time constant of 30 msec is about 0.75 of the initial voltage. 


To illustrate the graphical interpolation for two variables consider the above 
problem, except make the final rate 5x10° Rads/sec. Now we must interpolate for both 
the time constant and the final rediation rate. One way of doing this is to repeat 
the previous interpolation for several final rates and plot the voltage vs final 
rete for a time constant of 30 ms. The four curves in Figure 5 show the interpolation 
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for the time constant of 30 msec for final rates of 10", 10°, 10°, and 10! 


Reds/sec. Figure 6 then plote the four voltages for a 30 msec time constant 
ve the final rate. From this figure we see that the voltage left is about 0.52 
of the original voltage. 


Obviously this same technique can be used for any mamber of variables, and 
the order in which the interpolations are made is not too critical. The number 
of points used is strictly a mgtter of judgement. I used four here simply to il- 
lustrate that two might not be enough. Of course, there is no guarantee of accuracy 
in these interpolations, but in view of the expected capacitor-to-capacitor response 
variations, the accuracy should be quite sufficient for engineering use. 


Another problem which might arise in the use of these tables is how do we 
calculate the voltage left if the initial pulse has a different length than the 
25 nsec used for the initial pulse used in these tables? The answer is, very 
simply, that if the initial pulse length is short compared to all capacitor and 
circuit time constants, then the voltage at the end of the pulse is proportional 
to the duration of the pulse. A good engineering approximation to the voltage 
after any length of time can be made by plotting the case of interest from the table 
with the initial pulse length equal to 25 ns and then plotting initial drop for the 
other initial pulse length and making the voltage left at longer times assymtotic 
to the first plot. An easy interpolation for the initial pulse magnitude may be 
found from the fact that the voltage drop from the first pulse is proportional to 
the delta'th power of the initial radiation rate. 


3.2 Resistance Simulation Tables 


The purpose of these tables is to help the designer simulate the effects of 
the radiation pulse in his circuit without actually irradiating the circuit. Con- 
sider the switch in Figure 3 to be a relay and R, to be a fixed resistor. Then ve 
might ask the question: for a particular capacitor type, a particular radiation 


‘rate » & particular circuit time constant, and a particular time at which we are 


interested in the performance of the circuit, what value of resistance should ve 
use in the above simulation to give the same capacitor voltage at the time of in- 
terest that the radiation pulse gives? The resistance simulation tables are de- 
signed to help answer this question. Each table is calculated for a different time 


ana avpestissestintatesiat srnsiats ieee onleee 
SCH Et Be 


4 ea arate fa 
HaULSEETGTSUUNGGIEGHVEATE HERE 
a iil om 


Sean Este HEEGatTE quit 
RTH Hat ttt a Ur iD, 


I 2 ia We 
pees teais 


seeneiattitar HE i iat 13 
H seseed 
et EE tiie a | 
Hiatt Hy 3 


Hert Lae 
os 


SGuGcuam) wouccuan SSRSRTR SETS Mth Hoth 
SSSKHED (AOSERESRAGSUAOSTSGRRBABEE A AFRER 
am asabense couse 


meGuRA Setietel tecaeoes eet aunea 
a 


acacre 
MP? on 8 
ZENFET SI salem cS 


MEPaaEr Fri 2t 


° 


Bee 


suaNe Gusecesa veaasa sosa = 

rr tt et 

SHRRERUURRSASET ARE eeeeaee a 
seu sesnepeseseegeees 


as é Soo Gans 
tp eT : . x + * 
frgby peagiit pesthedtat ytestenaies 77 cee 
SS ESRES EESIOEEECIENEES oR a Br: -JESEZESSEE aati 4 3 “| 
b EERSTE SESGENRESESZEEUEOSES: SHEESH EAGLE le 
re En an rr 
SLES SEE HITE Hi 


wee ee a + pene eee — ween a =e ooo mee ees 


sEaRRRnaooeeed aes 


BHaE 
oe et 


‘Ginna Ha a ia 
3 ET 


, 
> 
} ; ‘ ‘ 
oe i 
f Hil 
- < dab 4 
' * , | 
4 -+ 
; 


ve — 


gausaes queeste 
ABESSURLES 
S ERASE USES EDSSSAR 


AK e ATT HH 

sete ESSUEEHVESLHEGEET a 
: BNGAATE FREI VOENIe<QHII09? 
CUNT nee 
fast oa 


*- 
—_ 


ie 


HL AT Hi 


imate MLLER LL 
a aia a8 |e 


ce 


5 iiaie py an Fewatauea tneeutvege ovation 
tL a tt iw alae , 
a init 
i it Dann HL dN 


oa 


| POLE UTP 


SRST FCM (HB ifih | 


a 
SSRCAAN Ht 


) 
-—= ee 


Me votes bad i we ‘ = t 
iat ra aan Lg 
TH Pay Sais aan iasb-t 


acacrca 


PsO Bers se 
ee 


¢ . 
32a! 


WEN AE? O CAete 
SEW TOCSY AMS 


K-= 


i) 


32BG8 Ssebancunne Gnawa wamueepause Sanna 


fgets 


Th gus pgeaasuear reset senne 
SSTEdREETACCSUOSETOREROROGSEQUEE 
sat cn 


an Sn 


fit att PSeEEOREEE Bt I 
Sania ae if thoi 
SUERUSER20SGR00 20008 0008 Try 


in Li f 


sasae it A isco 
anges Gunns svwneceuuguacsn q sean ark 
Petit. 
Gaasneet reece Reeeteentte Cal 
ee | ee Oe 
poo Htc Water 
pee HE tH Hf ii ‘ia ee ay Pa 
atte lt : 
HTT ee ae IH Lite 
aE AEE EE sa ee ETE 
cH pa ee Shon oan ee 
suaasnsesusr, 8 i SOA ee ae bn LEN PETE AN 
pesGusauedg cease coneas | BREN WS) DEY” 
SUC TRASS sicsstsrsstsensteasvesce at pthieebieh sbectotasl Serer 
SERSAURERE SSRSREAERY ARSENE SCY, a Ce PRL 
| Ben i ws Se 
% os Th io-0.) Rieti Hen Wkt lonad ; 
wetatwes 
EA |B it hy Oo A Ee 
iittala RE HME 


—_— 
44} 
pombe e ‘ 
oe t 
: i ; 
' in | 
' } ‘ ‘ 
ih tf 
—-+ 
’ ' ‘ 
Tit 
i ! i 


oN + ---+- + 

1. Seukena -_— > ~- \ 
Foe wMa Seseteeceeeeeeae PE 1 Ra 1 WR ME me 
HanELLSALELES Seeeena SERARSERES PRUE BEER Looe anil 
siensessreretrsee ptf mane siti is TS 

sreqataenuesastense ee PEON Wonks 

Banaaueaes i eo 

‘ 


ttt 
Seiiscie cn Ele HH iia 


t H 
BEESSES222c 2829293 EPEC FiSaTEE rv 
HES UAT 


oe ee oe ee ee ee oe ee ee 
Oe pe re ee ee ee ee ee ewe: 


after the start of the radiation pulse. Across the tops of the tables are 
various circuit time constants, corresponding to various R,'8 in Figure 3. 

The row headings are the time constants corresponding to R, in Pigure 3. The 
use of the tables is: WMirst, look up the capacitor voltage left at the time 
of interest in the capacitor tables. Then go to the resistance table for the 
time of interest and find the column corresponding to the circuit series time 
constant. Follow this column down until the voltage shown is about equal to 
the voltage found from the capacitor tables. The time constant at the left end 
of the row in which this number is found gives a value of resistance to use in 
a circuit simulation. 


As an example of this technique consider that we wish to simulate the dis- 
charge of a mylar capacitor ina circuit with a series time constant of 10 msec. 
The radiation environment is an initial rate of 1012 Rads/sec and the final rate 
is 10° Rads/sec. The time of interest is 100 msec. From the appropriate mylar 
capacitor table we find that the voltage left under these conditions is 0.852 of 
the initial voltage. Going to the resistance table calculated for 100 ms and 
following the 10 ms time constant column down, we find that an equivalent fixed 
resistance for use in a simulation circuit corresponding to Figure 3 would be one 
that gave a time constant of between 1 and 10 ms and closer to 10 ms. If the capaci- 
tor were 1 ufd, then the equivalent resistance would be between 1000 and 10,000 ohms, 
and closer to 10,000 ohms. 


4.0 ACCURACY OF TABLES 


The resistance simulation tables are exact to the number of significant 
figures given. It is believed that most of the numbers in the capacitance tables 
are accurate to three significant figures assuming absolute parameter accuracy. 
There are some erroneous points in these tables however. These errors arise be- 
cause the numerical integration used in the computer program may run into inaccura- 
cies under certain conditions. The integration routine used was a 201 point Simp- 
son's rule procedure. There were no iterative or checking loops since experimen- 
tation showed that 200 points were enough to make most points accurate and check- 
ing would add unconscioneably to the computer time. Errors in the table are fairly 


= 23 


obvious. They may arise whenever 1/200 of the time at which the number is cal- 
culated is on the same order of magnitude as, or greater than some circuit or 
capacitor time constant, but the capacitor has not yet reached a steady state. 
Those numbers satisfying the above conditions which are in error may generally 

be recognized because the rate of change of voltage in the table is faster than 
seems physically reasonable. Of course, variations in capacitors will give larger 
errors than any in the tables. Also, most of the parameters were measured at rates 
mach lower than the initial rates used here. Extrapolation of data to the higher 
rates may also result in errors. 


5-0 MATHEMATICAL DERIVATION 


Consider the circuit in Figure 3. Assume that at time t = 0 the switch is 
closed. Then, from Kirchoff's current law: 


i,2ig+ i, . (1) 


By the definition of a capacitor: 


| Ye 
4, ° Car» (2) 
and, from Ohm's law: 
V 
c 
i, = R, ’ (3) 
Vi-V 
Ss = oan eee r (4) 
0 
Thus: 
vV.-V Vv av 
a o 
re) 8 
If we set: 
R=,» (6a) 
CR, s tT.) (6b) 
Vv 
a wx, (6c) 
° 


cn ae ne ee en a EE Pn meet , 1 - 
wh « 


- 1h - % 


then equation 5 becomes: 


T 


tre tT) Fz - (7) 


Equation 7 is the basic differential equation which must be solved to get 
the transient response of the capacitor circuit. Assumptions listed in Section 
2.1 will be used. This means that any time variation in Equation 7 must be in 


the T. term. 


From the TREE Handbook, T. is a function of time and is given by: 


t 
A 
FP) oR P(e) ey | on [s- : J Menace’ 
dl 


' ‘ 
+ Kao on = J y Acerat” ’ (8) 
Tao 
@ 
vhere Ky 4, Kay? Tay? Kao» and Tap are parameters of the material used in the 
capacitor dielectric, and y(x) is the rediation rate in Rads/sec at time x. 
Equation 7 may be written to show its time dependence as: 


mtxje + f(t] eS. (9) 
° 


Thies is a linear, first-order differential equation whose solution may be written: 
t 


x(t) = ae op ( [2+ x] «| x 
6 


aa -f[R- *(t)] at) , (10) . 
silo 
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where G is an albitrary constant selected to make x = x, when t = 0; i.e. if 
we set: 


F(t) = exp ( st + r(+)| at) , (11) 
fe) 
then, at time t,: 
t 
“ alt,) = (gees) (xpr(0) +b | “r(eyat) . (22) 
) F(t.) ° ie : : 


Equations 8 and 1) may be solved to get F(t) in closed form. Two solutions 
are necessary, one for times during the initial pulse and the other for times dur- 
ing the final pulse. These solutions are incorporated in the computer program and 
are used to mmerically integrate Equation 12, which has no closed form solution. 


Equation 7 may be solved in closed form if T. is constant. This solution 
is, forx = 1 at t = 0: 


ae L exp |= h ney 
To) Ts T) T. (13) 
ene 
. 7. t. 


The resistance simulation tables were calculated by substituting various values of 
To Te» and t in Equation 13. 
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APPENDIX: THE COMPUTER PROGRAM 


The computer program used to calculate the capacitor tables was written 
in the version of ALOOL used by the Burroughs B5500 computer. A listing of the 
’ program is given in Figure Al, and a general flow chart is given in Figure A2. 
A description of the programming language may be found in Burroughs Corporation 
Report $5000-21012, "Extended ALGOL Reference Manual for the Burroughs 35000." 


There are a few points which should be made for the benefit of anyone who 
wishes to rewrite this program for another machine. 


1. A lot of computer time can be saved by minimizing the number of 
exponentiations required in the calculations. One way to help 
attain this goal is to make use of the fact that e™ = (e*)" as 
is done in several places in this program. 


2. If integration or differential equation solving techniques are used 
which adjust their step lengths until the errors are sufficiently 
small, then each table would probably require hours of computer time. 
This is because the step length will have to be less than the shortest 
time constant in the problem. The technique used here is to notice 
from Equations 12 and 8 that x approaches a constant at long times, 
to replace divergent values vith this constant, and to accept the 
few errors that occur. Computer time for one table was about 20 
sec, On the faster machines, say like the IBM 7094, it is estimated 
that one table would take about 10 sec., exclusive of I/O time. 


3. Care should be used in factoring cerms like 
g(at - dt) 


into 
ett et 


as it is easy to run into situations such as 


6 
0  ,-20 
* in these calculations, giving magnitudes which no computer can handle. 
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V@” 
BEGIN COMMENT TRANSIENT RADIATION CN CAPACITORS PROGRAMS 
INTEGER Wo 135 Ilo KKB 
COMMENT WN IS NUMBER OF INTERVALS FOR SIMPSONS RULE AND WUST BE EVEN 
II 18 THE NUMBER OF CAPACITOR TYPesS Jy 1S THE NUMBER OF INITIAL RAT 
ES KK 38 THE NUMBER OF FINAL RATES) 
READ(C Ns Ble JJo KK 
@EGIN REAL QO) | 
ARRAY ZAs 2000 & 2393 
REAL WAs wo wid 
ARRAY 222(0 8 8)) 
ARRAY RTIMs TIMCO & 23)5 OUTPs XO» TAURECCO &¢ 8)» INGR»e Eho EBs 
TTCO 8 No OLAs KPAs KAA» TAAs KOAsn TRACE § B3)5 INRAs TIMRACO 
@ dJJJe FINRACO & KKID 
ALPHA ARRAY ALLPCO & 23)5 ALFYCO 8 ff5 © & G@)o ALFRACO & JJ?) 
BOOLEAN BOLA, B80LB3 
INTEGER Io Jo Ko Lo Leo LLLD 
REAL KP» DEL» KAs TAs KO» TOs Sas SAo As Bo Co Oo AAs BBs CC» 
AAAs 880s CCCs OAs OB8s XT» VAs Vas 


FORMAT 
FOACKIGS "THIS TABLE IS THE CALCULATED TRANSIENT RADIATION RESPONSE*// 
R320"OF A %5586// 
N2Go°THE CAPACITOR PARAMETERS USED ARES" 
M2to KPa", £11.45" KOLO rEg ,Oo" KO20%sE11.,4/ 
KIGoPDEL TAS SE11,00% TD018%,E11,49" TO2E"5E11,4// 
NISso"°THE RADIATION PULSE £8"%rE19045" RADS/SEC FOR"DE11.45" SEC*/ 
N20s"FOLLOWED BY A CONSTANT RATE OF 90611040" RADS/SEC*S/ 
MISs"TAU IS THE TIME CONSTANT OF THE CAPACITOR CHARGING CIRCUIT.” J)» 
ee eae VALUES GIVEN ARE THE RATIOS OF THE CAPACITOR VOLTAGE AY TH 
N22o°TIME ZSNOICATED TO THE INETIAL CAPACITOR VOLTAGE" /)>, or 
FOCCK2Z2s"TAU (TIME CONSTANT)®/ 
MAD oP IME MKS 0910 SECM ONO] SECMON259100 MSPoKGr 910 MS DX» | 
PE MSF aK 907100 US%XG0%10 USM en XSo%1 USP /"TIME"), 
_FOOCAG) X309CF6.40K3)99 . 


wee" FILE RUIMCOD WETH Os Or 0965 20265 S@265 C955 2095s S095, Body 


2020, SOR, B93, 20P35 SOP3, 092s 20°92, SOP2s Fein 20°15 5m Ip e 
20 Se 103 

FILL ALLPCe) WITH * 0 My MUXKKKKMs © § US®s © 2 US", 

# § US®> * 10 US®s * 20 US*» * $0 US*s 7400 US"» "200 US">» 

"$00 US"» © 1 MS%> * 2 MS%s © GF MSM, * 10 M§%> * 20 HE» 
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FILL TAURECECe) WITH Oo to Io 105 1005 10005 OA» OSs 063 © 
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TTCO) ¢ OD 
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REAODCFOR f © 1 STEP 1 UATIL II NOCFOR yo 0 STEP 1 UNTIL &@ CO 
ALFYC IT» J)o DLAC TL )o KPACT I)» KAACT I» VAACT)» KBACI)» TAACII)) 
READCFOR IT ¢ 1 STEP $ UNTIL Jd DOLCINAALCT I» TIMRACT I» ALFRACTII)): 


REANCFOR 3 ¢ 1 STEP 1 UATEL KK 00 FINRACII)) 
FOR L ¢ 1 STEP 1 UNTIL IF 00 
BEGIN KP © KPACL 9) 
DEL ¢ OLACLID 
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‘ KA @ KAACLIS 
ot Ka ROALD 
ded TA © TAACLYD 
} Te ¢ TOACLIS 

BOLA ¢ KA = OF 


BOL8 ¢ KB s 03 
$A © KA & TAS 
$8 ¢ KB x TBs 
, FOR LL ¢ &§ STEP 4 UNTIL JJ 00 
BEGIN TIME) © RYIME4) © TEMRACLLID 
ALLPC1) © ALFRACLL ID 
ZAC1) © IF BOLA THEN 2.0 ELSE ExPCe Tree 
2001) ¢ IF BOLO THEN 1.0 ELSE EXPCoe Tine 
FOR I ¢ 2 STEP 4 UNTIL 23 00 
BEGIN TIME 3) © RIIMET) © TIMCL)S 
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BEGIN OB © FINRACLLLI © OEL3 
AAA © OB RCKP o SA © $899 
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© TEIMC1) 4 1809993 
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FOR 1 © 4 STEP § UNTEL 23 00 
BEGIN WY © TEMCEI3 
WA @ ZACI}) 
we ¢ 20019) 
TTC1) © RT oOo WT 7 OND 
VA © EACI) © IF BCLA THEN 1.0 ELSE EXPC® RT / TADS 
VO © EOC!) « IF OCLA THEN 1.0 ELSE EXPCe RT / TO)3 
FOR y © 2 STEP ¢ UNTIL WN 00 
BEGIN EalJ) © EACY © 139 m VAD 
Eat’) © €fY © 1) = VO) 
TTC J} © THCY © 13 4% KTS 
ENO FILLING OF TIVE ARRAYS) 
FOR y © 0 STEP 1 UNTIL NW OO INORE YD © EXPCR nCEAC YU) 
@ WA) © C BCEBC I) © WODD) 


Pe —« ° . wiv en “s 
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FOR yy © O STEP 1 UNTIL 4 00 
OECIn G@ © CUTPCU)) 
O  TAUREC( JI) 
' Eatn) ¢ 1.03 
Va © EACKR © 1) © EXPCe KT CO © ADDS 
FOR K © A © 2 STEP © ¢ UNTIL O 00 EACK) © EACK + 
1) = VAS 
VA © EACO) x IAGREIO) » EACH) @ INGACNYS 
VB ¢ O83 
FOR K ¢ 2 STEP 2 UNTIL N © 1 00 VB © V@ 4 INGREK 
) « €ACKID if 
VA © 4.0 = VB ¢ VAS 
ad ve ¢ O39 
FOR K «© 2 STEP 2 UNTIL N © 2 00 VO © VO + INGRIK 
) = EACK)S 


Va @¢ X¥ =» DO e(200 » VO @ VAD 4 9) 
OuTPtJ) © XOrJ) =» EXPce wt (0 6 AD) > B HCE @ hA 
> C€ Cd 7 WB) © VAD 
If 1.0 < QUTPCJ) THEN OCUTPC I) © Z22CJ) ELSE IF 2 
< 1 AND CUTPEJ) > ZZZ2CN) ANDO OC S$ ZZ2Z20U) THEN 
OuTPeJd) © ZZ2ZCJ)) 
ENO FIX TIME VARY Taus 
WRITECFHDOs ALLPCI)» FOR y © O STEP 1 UNTIL 8 OO 
GUTPe))) 
if 1 8 ( THEN 
BEGIN A © AAAS 
8 ¢ BBA) 
C ¢ CCC 
FOR y ¢ O STEP 3 UNTIL 6 00 XOCy) © OUTPC YS 
END COEFFICIENT CHANGED 
; END STEP THROUGH TIMES ‘ 
: WRITEG PAGED) 
ENO FINAL RATES) 
ENO INETSAL RATES) 
END CAPACITORS) 
ENO BNNER BLOCKS 
Eno. 


oe e@e: -m we o; wee CS ee} 2 e . oo 
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Read: n = # intervals for Simpsons 
rule integration 

ii = $ capacitor types 
-Jj = $ initial radiation rates 
kk = # final radiation rates 
capacitor parameters (11 sete) 
initial rate info. (jj sets) 
final rete info. (kk sets) 


Initialize variables. Calculate all 
numbers which will not change with 


cap. type, int. rate, final rate, 
time constant, or time. 


Select next capacitor type. 


Calculate all numbers which do nc 
change with int. rate, final rate 
time constant or time. 


Select next initial rate. 


Calculate all numbers which do not 
change with final rate, time 


constant, or time. 
Select next final rate. 


Calculate all numbers which do not 
change with time constant or time. 


Write: Table headings; row for 
time Ss 0. 
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me lire aS 


_—— 


Select next time. 


La. numbers, inclu 
parte of integrand, which do 
change with time constant. 


Select next time constant. 


Pinish Simpson's rule integration; 
Calculate output number. 


Have all time constants been used? 
yes 


Write: Output row for time used 
in calculations. 


be. Was calculation for initial rate 
only? 


nO_— Have all times been used? 
yes 


Skip to top of next printer page 
for new table. 


to 
3 


go 
Page 
ll 
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ore “ Bave all final rates 
: veen used? 
yes 
he Ee se Bave all initial rates 
<a deen used? 


yes 


~ S BO Have all capacitor types 
been used? : 
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